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ABSTRACT 

Using a large sample of local galaxies (144,940) with -17.5 < Mr <-22 and 0.025 < z < 0.107, 
selected from Sloan Digital Sky Survey Data Release 5, we compare AGN host galaxies with non-AGN 
galaxies at matched luminosity, velocity dispersion, color, color gradient, or concentration index, to 
investigate how AGN activity is related with galaxy properties. The AGN sample is composed of Type 
II AGNs identified with flux ratios of narrow-emission lines with signal-to-noise ratio > 6. We find 
that the fraction of galaxies hosting an AGN(/agn) depends strongly on morphology together with 
color, and very weakly on luminosity or stellar velocity dispersion of host galaxies. In particular, /agn 
of early-type galaxies is almost independent of luminosity nor velocity dispersion when color is fixed. 
The host galaxy color preferred by AGNs is u — r « 2.0 for early-type hosts and u — r = 2.0 ~ 2.4 
for late- type hosts. This trend suggests that AGNs are dominantly hosted by intermediate-mass late- 
type galaxies because early-type galaxies with u — r k, 2.0 are very rare. We also investigate how 
the accretion power varies with galaxy properties. While the Eddington ratio ([OIII] line luminosity 
normalized by black hole mass) ranges over three orders of magnitude for both morphological types, 
late-type galaxies are the dominant hosts over all AGN power. Among late-type galaxies, bluer color 
galaxies host higher power AGNs. These results are consistent with a scenario that more massive and 
redder galaxies are harder to host AGNs since these galaxies already consumed gas at the center or 
do not have sufficient gas supply to feed the black hole. In contrast, intermediate-mass, intermediate- 
color, and more concentrated late-type galaxies are more likely to host AGNs, implying that perhaps 
some fraction of low-mass, blue, and less concentrated late-type galaxies may not host massive black 
holes or may host very low-power AGNs. 

Subject headings: galaxies: active — galaxies: evolution — quasars: general ~ methods: surveys 
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1. INTRODUCTION 

Since the discovery of active galactic nuclei (AGNs) 
at the centers of nearby galaxies (Seyfert 1943), and at 
cosmological distance as highly energetic optical and ra- 
dio sources (Schmit 1963), it has been a long-standing 
issue why some galaxies host AGN while a large frac- 
tion of galaxies do not (e.g. Martini 2004; Jogee 2006). 
Are there fundamental differences between non-AGN and 
AGN host galaxies or are AGN host galaxies a random 
subset of normal galaxies? Spatially resolved kinemat- 
ics of very nearby galaxies revealed that massive bulge- 
dominated galaxies host dormant supermassive black 
holes at their center (Kormendy & Gebhardt 2001; see 
Ferrarese & Ford 2005 for a recent review) , implying the 
ubiquity of black holes in massive galaxies. Hence, trig- 
gering AGN activity is related to the fuel supply mecha- 
nism, by which dormant black holes become active. 

In turn, as a feedback mechanism AGN activity can 
play an important role in suppressing star formation and 
gas cooling in galaxy or galaxy cluster scales as demon- 
strated in various semi-analytic models (Croton et al. 
2006; Robertson et al. 2006; Bower et al. 2006). Under- 
standing how AGN activity is related with galaxy prop- 
erties and large scale environment may shed light on the 
origin of the black hole - galaxy connection, as observed 
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in the present-day universe (e.g. Ferrarese & Merritt 
2000; Gebhardt et al. 2000), and its evolution (Woo et 
al. 2006, 2008; Peng et al. 2006; Younger et al. 2008). 

Some galaxies with particular properties, indicative of 
an efhcient fuel supply mechanism, may present more 
AGN activity than other galaxies. In a standard galaxy 
evolution picture, AGN activity is naturally expected 
through galaxy merging and interactions, which can pro- 
vide gas to the center and feed the nuclear black holes 
(e.g. Sanders et al. 1988; Hopkins et al. 2006). Host 
galaxies of AGNs are expected to be different from non- 
AGN galaxies in this scenario, and this picture may be 
relevant to explain highly energetic AGNs at z ^^2-3 near 
the peak of quasar activity and star formation (Hasinger 
et al. 2005). Internal processes such as bar driven gas 
inflow (e.g. Combes 2003; see Kormendy & Kennicutt 
2004 for a recent review) , turbulence in inter-stellar mat- 
ter (e.g. Wada 2004), and stellar wind (e.g. Ciotti & 
Ostriker 2007) are also considered as a gas supply mech- 
anism, responsible for triggering AGN activity even with- 
out invoking galaxy mergers. 

However, it has not been clear whether AGN host 
galaxies have distinct properties compared to non-AGN 
galaxies. Various studies indicated that AGN host galax- 
ies do not show excess of bar presence or local density, 
compared to non-AGN galaxies (e.g. Combes 2003; cf. 
Maia et al. 2003). HST imaging studies of quasar host 
galaxies at low redshift showed that host galaxies have 
structures and stellar populations very similar to those 
of non-AGN galaxies (e.g. Bahcall et al. 1997; McLure 
et al. 2000; Dunlop et al. 2003). Photometric and 
spectroscopic studies of radio-loud AGNs showed that 
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their host galaxies he on the same fundamental plane 
as quiescent galaxies (Bettoni et al. 2001; Barth et al. 
2003; Woo et al. 2004, 2005). In contrast, the presence 
of cold molecular gas and young stellar population in 
quasar host galaxies suggested that these host galaxies 
are not typical early-type galaxies (Scoville et al. 2003; 
Tadhunter et al. 2005). A recent study with deep HST 
images revealed that some host galaxies have weak merg- 
ing signatures, indicating recent merging or interactions 
although majority of stars are relaxed as in the case of 
quiescent galaxies (Bennert et al. 2008). 

The Sloan Digital Sky Survey (SDSS) provides a large 
sample of non-AGN and AGN host galaxies, enabling to 
overcome the difficulties caused by the relatively small 
sample size in the previous studies (e.g. Ho et al. 1997; 
Hunt & Malkan 2004). Since the SDSS data became 
available, various effects on AGN activity have been in- 
vestigated, including host galaxy properties (Kauffmann 
et al. 2003; Rich et al. 2005; Schawinski et al. 2007), 
local density, and large scale environment (Kauffmann 
et al. 2004; Constantin & Vogeley 2006; Constantin et 
al. 2008). By studying a large sample of emission line 
galaxies, Kauffmann et al. (2003) reported that early- 
type galaxies with AGN have bluer color than those with- 
out AGN, indicating that AGN host galaxies are inter- 
mediate objects between the red sequence and the blue 
cloud (see also, Schawinski et al. 2007). Similarly, it 
was reported that early-type host galaxies of LINERS 
have younger stellar population than quiescent early-type 
galaxies (Graves et al. 2007). 

Since early-type and late-type galaxies are in differ- 
ent evolutionary stages as demonstrated in the color- 
magnitude diagram (e.g. Bell et al. 2004), and the 
galaxy mass is the main driver of galaxy evolution as 
observed in galaxy downsizing (Cowie et al. 1996; Treu 
et al. 2005; De Lucia et al 2006), it is crucial to com- 
pare AGN and non-AGN galaxies with matched galaxy 
morphology and mass scale. Then various effects due to 
galaxy evolution can be cancelled out. However, a de- 
tailed study with a large sample of matched host galaxy 
morphology and mass has not been reported. For a mag- 
nitude limited sample of nearby ~ 500 galaxies with de- 
tailed morphology information. Ho et al. (1997; see also 
Ho 2008 for a review) investigated the AGN fraction of 
each morphological type, reporting '--^50% AGN fraction 
in early-type (E-Sab) galaxies. The high AGN fraction 
in this study is due to the fact that high spatial resolution 
(< 200 pc) was obtained and, hence, faint nuclear emis- 
sion lines were detected without strong contamination of 
host galaxies. With a much larger sample of emission line 
galaxies from SDSS DRl, Kauffmann et al. (2003) re- 
ported that AGNs are preferentially found in more mas- 
sive galaxies and that host galaxies of high [OIII] lu- 
minosity AGNs have younger stellar age than non-AGN 
galaxies. To achieve more detailed understanding on trig- 
gering AGN activity, it is necessary to have morpholog- 
ical classification and matched host galaxy properties in 
comparison, and proper normalization of AGN luminos- 
ity by black hole mass. 

Using a large volume-limited sample with detailed 
galaxy properties, we initiated an intensive project to in- 
vestigate environmental effects on AGN activity in three 
different scales: host galaxy, small scale environment, 
and large scale environment. In the first paper of this 
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Fig. 1. — The volume-limited sample. The gray and blue dots 
in the whole area enclosed by rectangular boxes represent the total 
sample and the AGN sample, respectively. Each subsample with 
various volume is defined with rectangular boxes. The green curve, 
with Petrosian r-magnitude (mpet = 17.6) corresponds to the flux 
limit of the spectroscopic survey. 



series, we explore how AGN activity is related with host 
galaxy properties using a large sample of SDSS galaxies. 
Small and large scale environmental effects on AGN ac- 
tivity will be addressed in the next paper (Y.-Y. Choi et 
al. 2009, in preparation). We discuss our sample selec- 
tion criteria in § 2 and sample properties in § 3. In § 4, we 
explore how AGN fraction changes depending on the host 
galaxy properties. In § 5, we present how AGN power 
is related with each host galaxy property. Summary and 
discussion will be given in § 6. Throughout this paper, 
we adopt a flat ACDM cosmology with fl^ = 0.27, and 
flA = 0.73. Magnitudes are given in the AB system. 

2. DATA 

2.1. Sample Selection 

We first selected galaxies with the r-band absolute 
magnitude Mr < —20-1- 5logh (hereafter we drop the 
+51og/i term in the absolute magnitude) and 0.025 < z < 
0.107 from the large-scale structure sample, DR4plus, us- 
ing the New York University Value-Added Galaxy Cat- 
alogue (Blanton et al. 2005), which is a galaxy subset 
(the Main galaxy sample) of the SDSS Data Release 5 
(Adelman-McCarthy et al. 2007). Readers are referred 
to Choi et al. (2007) for a detailed sample description. 
The rest-frame absolute magnitudes of individual galax- 
ies are computed in the r-band, using Galactic redden- 
ing correction (Schlegel et al. 1998) and if-corrections 
as described by Blanton et al. (2003). The mean evo- 
lution correction given by Tegmark et al. (2004) is also 
applied. Our survey region, with angular selection func- 
tion of more than 0.5, covers 4464 deg^ (see the sur- 
vey boundaries in Figure 1 of Park et al. 2007). To 
increase statistics toward low luminosity and low-mass 
scales, we added five volume-limited samples selected 
from the same Catalogue, with fainter magnitude ranges 
down to Mr ~ —17.5. The definitions of all 6 volume- 
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TABLE 1 
Volume-limited Samples 



Name 


Magnitude 




Redshift 




Distance 


Number 


d 


sub6 


-20.0 > Mr 




0.025 <z < 0.10713 


74.f 


) < R < 314.0 


80,478 


5.56 


sub5 


-19.5 > Mr > - 


20.0 


0.025 < 2 < 0.08588 


74.C 


) < i? < 252.9 


31,048 


4.58 


sub4 


-19.0 > Mr > - 


19.5 


0.025 < z < 0.06869 


74.C 


) < < 203.0 


16,772 


4.18 


subs 


-18.5 > Mr > - 


19.0 


0.025 < z < 0.05485 


74.f 


) < R < 162.6 


9,153 


3.78 


sub2 


-18.0 > Mr > - 


18.5 


0.025 <z < 0.04374 


74.f 


) < i? < 129.9 


4,932 


3.41 


subl 


-17.5 > Mr > - 


18.0 


0.025 < z < 0.03484 


74.f 


X i? < 103.7 


2,557 


3.00 



Note. — Col. (1): name of subsamples. (2): absolute magnitude range in the r-band. (3): redshift range. (4): distance range in units 
of /i~^Mpc. (5): the number of galaxies. (6): mean separation of galaxies in units of h~^Mpc. 
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Fig. 2. — Distribution of galaxies in the emission line flux ratios 
plane. The AGN sample is defined by the solid curve (Kewley et 
al. 2006). To minimize contamination of star- forming galaxies, 
composite objects below the solid line are excluded. Blue dots 
represent high-power AGNs {i[oiii]/.'^^BH> 1 in Lq/MQ) while 
red dots represent low-power AGNs. 

limited subsamples are summarized in Table 1. 

The total sample of 144,940 galaxies combined with 
all subsamples is presented with gray dots in rectangular 
boxes in Figure [TJ The difference in the volumes of each 
subsample has been taken into account throughout our 
analysis. 

2.2. AGN Selection 

Since we will compare AGN and non-AGN galaxies 
with matched luminosity, velocity dispersion, color, color 
gradient, and concentration, we selected an AGN sam- 
ple using Type II AGNs, for which host galaxy prop- 
erties can be determined in the same way as in non- 
AGN galaxies, without suffering further complication 
due to the presence of a bright nuclear. Type II AGNs 
can be separated from star forming galaxies based on 
the flux ratios of Balmer and ionization lines (Baldwin, 
Phillips & Tcrlcvich, 1981; Vcilleux & Osterbrock 1987). 
Figure [2] presents the line ratios of emission- line galax- 
ies, for which all four emission lines ([OIIIJA5007, H/3, 
[NIIJA6583, and Ha) are detected with signal-to-noise ra- 



tio (S/N) > 6. Out of the 144,940 galaxies, 46,520 galax- 
ies (32%) satisfy the S/N criterion. Following Vcilleux & 
Osterbrock (1987) and other studies, we used three clas- 
sification schemes including low-ionization species. In 
practice, we selected AGNs using a conservative AGN 
definition from Kewley et al. (2006); 

0.61/(log([NII]/Ha) -0.47) -f 1.19 < log([OIII] /H/3), 

0.72/(log([SII]/Ha) - 0.32) -1- 1.30 < log([OIII] /H/3), 

0.73/(log([OI]/Ha) -t- 0.59) -f 1.33 < log([OIII] /H/3). 

We also used Kauffmann et al. (2003) demarcation line 
(dashed line in Fig. [2]) to identify composite objects, 
which contain AGN as well as extended HII regions; 

0.61/(log([NII]/Ha) - 0.05) -f 1.30 > log([OIII] /H/3), 

0.72/(log([SII]/Ha) - 0.32) + 1.30 > log([OIII] /H/3), 

0.73/(log([OI]/Ha) + 0.59) + 1.33 > log([0III]/H/3). 

We excluded potential Type I AGNs. SDSS spec- 
troscopic pipeline does not classify broad-line AGNs as 
'galaxies'. However, some narrow Type I AGNs can 
still be present among galaxies. We assume galaxies 
with a Ha emission line width larger than ~500 km s~^ 
(FWHM) as Type I AGNs. Only 412 objects were found 
in our sample and excluded. 

Combining pure AGNs and composite objects, we de- 
fined an AGN sample, which is composed of 11,521 ob- 
jects. The fraction of each class is shown in detail in 
Table d The AGN sample is 7.9% of the total galaxy 
sample, however we note that this fraction depends on 
the S/N criterion on the emission lines. For example, 
when we used S/N > 10, the number of AGNs drops to 
6,088. Thus, the AGN fraction presented in this paper 
should be taken as a lower limit. 

The emission line fluxes and stellar velocity disper- 
sion were measured using an automated spectroscopic 
pipeline, specBS version 5 (D. Schlegel et al. 2009, in 
preparation). The basic technique used in this pipeline 
is described by Glazebrook et al. (1998) and Bromley 
et al. (1998). It determines radial velocity, redshift, 
and spectral classifications by fitting the spectra with 
SDSS-derived stellar templates as well as stellar tem- 
plates drawn from the high-resolution spectra from the 
ELODIE survey (Prugniel & Soubiran 2001). The li- 
brary of stellar templates is constructed by applying the 
principle component analysis technique to a sample of 
pure absorption-line galaxies. Note that these emission 
line measurements are slightly different from those given 
by Tremonti et al. (2004) because of the difference be- 
tween algorithms in subtracting stellar continuum be- 
fore carrying out emission line fits. However, only H/3 
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line fluxes are somewhat affected by this steUar contin- 
uum correction, with ^ 0.13 ± 0.10 dex difference in log 
([OIIIJ/H/3) compared to Tremonti et al. (the differ- 
ence in log([NII]/Ha) is only 0.02 ± 0.02 dex). Thus, 
a small fraction of composite objects in our sample was 
classified as star-forming galaxies in Kauffmann et al. 
(2003). However, our AGN sample with S/N > 6 is still 
more conservatively defined than that of Kauffmann et 
al.(2003) with S/N > 3. 

We adopt the [OIII] line luminosity (L[oiii]) sls an ac- 
cretion rate indicator, after correcting for dust extinction 
by using the relation (Bassani et al. 1999); 

i[oiii] =i[OIII],obs[(Ha/H/3)obs/(Ha/H/3)o]'•^^ (1) 

where an intrinsic Balmer decrement (Ha/H/3)o = 3.1 is 
adopted (Kewley et al. 2006). Black hole mass (Mbh) 
was estimated from stellar velocity dispersion assuming 
the Mbh-c relation (Tremaine et al. 2002); 

logMBH/Af© = 8.13 + 4.02 x log(CT/200 km s"^), (2) 

where a simple aperture correction to the stellar velocity 
dispersion has been applied (Bernardi et al. 2003). Then, 
the Eddington ratio is quantified by £[oiii] normalized by 
Mbh, as an indicator of AGN power. 

3. SAMPLE PROPERTIES 

To investigate how AGN activity is related with each 
host galaxy property, it is crucial to match other galaxy 
properties in comparing AGN and non-AGN galaxies. 
Morphology and mass are probably the two most fun- 
damental parameters of galaxies, presenting the current 
evolutionary status. By fixing morphology and mass, we 
can avoid confusions resulting from the correlations of 
other physical parameters with mass and morphology. 
In this section, we discuss our morphology classification 
scheme (§ 13. ip . the luminosity and velocity dispersion 
distributions 13. 2|) . and the Eddington ratio distribu- 
tions of the AGN sample (ij 13. 3p . 

3.1. Morphology Classification 

Morphology classification is not straight forward espe- 
cially for faint and seeing-limited distant galaxies. To 
apply to a large sample of SDSS galaxies, we adopted 
from Park & Choi (2005) a classification tool in the color 
(u— r) versus color gradient (A{g^i)) plane, where early- 
type and late-type galaxies are well separated. The g — i 
color gradient is defined as the color difference between 
the annulus with 0.5i?Pct < R < Rpct and the region 
with R < 0.5i?pet, where i?pct is the i-band Petrosian 
radius. The concentration index (cjn) was also used as 
an auxiliary parameter, which is defined as R^q/Rqq, 
where i?5o and Rgo are the radii from the center of a 
galaxy containing 50% and 90% of the Petrosian flux in 
the i-band, respectively. The boundaries between the 
two types in the three-dimensional parameter space is 
determined in such a way that the classification best 
reproduces the visual morphology classification: ellipti- 
cals/lenticular (early- type) and spirals/irregulars (late- 

type). 

Using this classification tool, we divide our sample into 
two morphological groups. Note that our morphology 
classification is mainly based on color and color gradient, 
however, the reliability of this scheme is ^^90% based 
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Fig. 3. — Segregation between early- type (left) and late- type 
galaxies (right) in the u — r color versus inverse-concentration index 
(upper) and u — r color versus and g — i color gradient (lower). 
Gray points are non-AGN galaxies, and blue (red) points are host 
galaxies of high(low)-power AGNs. [See the electronic edition of 
the Journal for a color version of this figure] 

on the test with training samples, for which color-based 
morphology was compared with morphology determined 
by eye (see Park & Choi 2005 for details). For practical 
reasons, we will call these two groups as early- type and 
late-type galaxies. 

Figure [3] shows our sample galaxies in the u — r color 
versus Cin plane (upper panels) and in the u — r versus 
A(5 — i) plane (lower panels). Red (blue) dots represent 
low (high) power AGNs as in Figure while gray dots 
denote non-AGN galaxies. Note that the sharp bound- 
aries shown in the u — r versus A(g — i) plane are due to 
the morphology classification criteria and that the scat- 
ter across the classification boundary is caused by the 
concentration index constraint (Park & Choi 2005). 

Our morphology classification yields a significant frac- 
tions of early types with bluer color and positive color 
gradient (bluer inner part). Among early-type galaxies, 
54% of AGN host galaxies is bluer than u — r = 2.4, 
while only 6% of non-AGN galaxies have u — r < 2.4. In 
the case of late- type galaxies, 17% of AGN hosts are red- 
der than u — r = 2.4, similar to that (15%) of non-AGN 
galaxies. These results imply that galaxies with inter- 
mediate color have higher fraction of AGN (see § 14.21 for 
more discussion). 

Table [2] presents AGN fractions for each morphology 
group. In the total sample, 7.9% of galaxies host AGN, 
of which 83% are late- type galaxies. 27% of AGNs has 
high accretion rate, L[oiii]/AfBH>l (Lq/Mq) and 84% 
of these high-power AGNs are hosted by late-type galax- 
ies. When we further exclude composite objects, 1.8% of 
galaxies host AGN, of which 72% are late-type galaxies. 

3.2. Luminosity and velocity dispersion 

Figure |4] presents the relation between luminosity and 
stellar velocity dispersion (a) for non-AGN and AGN 
host galaxies. For the stellar velocity dispersion, only 
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TABLE 2 

Sample Statistics 



% (Number) 


All 


Early 


Late 


All galaxies 


100.0 


41.3% (59,799) 


58.7% (85,141) 


pure AGN 


1.8 


27.5% (718) 


72.4% (1,887) 


composite 


6.1 


13.3% (1,191) 


86.7% (7,722) 


total AGN 


7.9 


16.6% (1,909) 


83.4% (9,609) 


Non-AGN 


73.5 


43.2 (45,981) 


56.8 (60,516) 



Note. — Col. (1): class. (2): fraction in total. (3,4): early- 
type galaxy fraction and number in each AGN class. (5,6): late- 
type galaxy fraction and number in each AGN class. AGNs are 
classified based on the ratios of emission lines detected with S/N > 
6. Note that non-AGN includes all galaxies excluding pure AGNs 
and composite objects defined by emission line ratio without S/N 
cut. 
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Fig. 4. — Luminosity versus stellar velocity dispersion of early- 
type (upper) and late- type galaxies (lower). Gray dots represent 
non-AGN galaxies while blue (red) dots represent high (ioi«)-power 
AGNs. Volume-limited subsamples are displayed with a magnitude 
gap for clarification. [See the electronic edition of the Journal for 
a color version of this figure.] 



Fig. 5. — [OIII] luminosity versus Mbh for AGNs hosted by 
early-type (left) and late-type galaxies (right). Dashed lines from 
upper left to lower right indicate ~100%, ~10%, and ~1% of the 
Eddington limit, respectively, assuming bolometric correction 3500 
for i[oiii] (with an uncertainty of ^0.4 dex; Heckman et al. 2004). 
We define high-power AGNs with i[oiii]/A^BH> 1 (Lq/Mq), cor- 
responding to >10% of the Eddington limit. Upper panels show 
relative fraction of AGNs as a function of Mbh- [See the electronic 
edition of the Journal for a color version of this figure] 

spectra with mean S/N per spectral pixel greater than 
10 were used. A simple aperture correction due to the 
finite size of optical fiber has been applied to the velocity 
dispersion (Bernardi et al. 2003). 

It is already noticeable that the distribution of AGN 
host galaxies is different from that of non-AGN galaxies 
in the sense that for given velocity dispersion, AGN host 
galaxies have higher luminosity than non-AGN galaxies, 
in particular for early-type galaxies. We will discuss in 
detail the relation between AGN activity and luminosity 
or velocity dispersion in § 14.11 We note that at the same 
luminosity, late-type galaxies can have a systematically 
wide distribution in mass. This may be due to the var- 
ious levels of young stellar population fraction. Thus, 
direct comparison between two morphological types is 
not trivial. 

3.3. Eddington ratio distribution 

Figure [5] presents distributions of AGNs in the i^[oiii] 
versus Mbh plane. Since Mbh correlates with the stel- 
lar velocity dispersion of the spheroidal component, the 
derived Mbh of late-type galaxies suffers systematic un- 
certainties. In practice, however, we cannot separate the 
bulge stellar velocity dispersion for our sample of SDSS 
galaxies, and take the total velocity dispersion to es- 
timate Mbh- Low Mbh values should be taken more 
cautiously. First, Mbh estimates lower than IQ^-^Mq 
are not reliable since velocity dispersion measurements 
lower than the instrumental resolution ('^ 70 km s~^) 
are not robust. Second, there is a deficiency of objects 
with .L[oni] ^ lO^i©, presumably caused by the detec- 
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Fig. 6. — The type-specific luminosity (panel a) and velocity 
dispersion [panel b) functions of non-AGN {dashed lines) and AGN 
host galaxies {solid lines), where the bin size of each plot is 0.5 
mag and 30 km s~^ for panels (a) and (b), respectively. [See the 
electronic edition of the Journal for a color version of this figure.] 

tion Hmit of the [OIII] hne since weak [OIII] hnes can 
be easily diluted by the host galaxy light. Therefore, 
we will focus on AGNs with Mbh > lO^ "^ M© (i.e. a 
^ 90 km s^^) in the following analysis when host galaxy 
properties are related with AGN power. 

We take further caution due to the fact that disk galax- 
ies can have shock-oriented extended emission line re- 
gions, which could mimic AGNs properties (Dopita & 
Sutherland 1995; sec discussion by Ho 2008). The po- 
tential host galaxy contamination is unavoidable in the 
SDSS spectroscopy with a 3" diameter aperture, corre- 
sponding to several kpc for typical galaxies in the sample. 
However, if the [OIIIJ/H/3 ratio of extended emission line 
region is different from the nuclear region, the contami- 
nation may not be severe (see Fehmers et al. 1994). 

We split our AGNs into high and low-power subsam- 
ples at i[oiii]A^BH= 1 (Lq/Mq), which corresponds to 
~ 5 — 10% of the Eddington limit (top dashed line) de- 
pending on the bolomctric correction (Hcckman et al. 
2004). Figure [5] shows that i[oiii] varies roughly from 
the Eddington limit down to ~ 10~^ of the Eddington 
luminosity (bottom dashed line). However, at the lower 
Mbh(< lO^'^Af©), the Eddington ratio range becomes 
narrower due to the incompleteness. Thus, we will use 
the high-power subsample (i[oiii]/-^BH> 1) to test any 

[OIII] "related selection effect in the following analysis. 

4. EFFECTS ON THE AGN FRACTION 

In this section, we examine how the fraction of galaxies 
hosting an AGN is related with host galaxy properties. 
We define the AGN fraction (/agn) as the ratio between 
the number of AGN hosts and that of all galaxies at fixed 
galaxy properties. We will investigate how /agn varies as 
a function of different galaxy properties: luminosity and 
velocity dispersion (§ 14. ip . color (§ 14. 2p . color gradient 
(§1131), and light concentration (? l4.4p . 

4.1. Luminosity and Velocity dispersion 

Figure [5] presents the number density distributions 
of non-AGN (dashed line) and AGN host (solid line) 
galaxies as a function of luminosity (M^) and veloc- 
ity dispersion. Morphology dependence of the distri- 
butions is examined by dividing the sample into early- 



type (circles) and late-type galaxies (stars). The shape 
of luminosity function and velocity dispersion func- 
tion of AGN host galaxies are qualitatively different 
from those of non-AGN galaxies. From the ratios 
(piE/L AGN)/(j){all AGN), one can see that AGNs re- 
side preferentially in luminous {Mr < —19) or intermedi- 
ate velocity dispersion (100 < a < 200 km s~^) late-type 
galaxies, although at higher velocity dispersions (> 200 
km s~^) early- type galaxies host comparable number of 
AGNs compared to late-type galaxies. Note that since 
the emission line strength roughly scales with black hole 
mass (and galaxy mass), detecting AGNs based on the 
emission line diagnostic is much more difficult at low- 
mass scales, which is reflected in the decreasing num- 
ber density of AGN host galaxies at lower luminosity 
(Mj. > —19) or lower velocity dispersion (cr < 100 km 
s-i). 

To understand the coupled dependence of /agn on Mr 
and a, we inspect how the /agn varies in the luminosity 
versus velocity dispersion parameter plane in Figure [71 
We present galaxies with cr > 90 km s~^ to reduce sys- 
tematic errors induced by the instrumental resolution of 
SDSS spectroscopy. 

On top of the distributions of the AGN host galaxies 
(blue dots) and all galaxies in the sample (gray dots) su- 
perposed are the contours of the constant AGN fractions 
in colored solid lines, while the dotted contours represent 
the number density of all galaxies in the sample. The 
smooth distributions of AGN fraction are obtained from 
the ratio of the sum of the weighted number density of 
AGN host to the sum of the weighted number density of 
all galaxies within the smoothing kernel at each bin in the 
parameter space. A fixed-size spline-kernel is used to cal- 
culate the weights. Contours are limited to regions with 
statistical significance above la for the fraction. The pa- 
rameter bin sizes for making the contours are denoted in 
each figure caption. A median uncertainty of the /agn 
is 5%, 9%, and 5% for all, early, and late-type galaxies, 
respectively. The /agn can be interpreted as the prob- 
ability for a randomly chosen galaxy to host AGN. As 
shown in Figure [71 the /agn is less than ^25% over all 
Mr or a ranges, regardless of morphological type. How- 
ever, we note that /agn is a lower limit since AGNs with 
weak emission lines (S/N < 6) are not included in our 
AGN sample. 

For the combined sample of early and late-type galax- 
ies in Figure [71 (left panel) , the /agn increases as host 
galaxy luminosity increases, suggesting that more lumi- 
nous galaxies are more likely to host AGN. In contrast, 
for given host galaxy luminosity, /agn typically increases 
with decreasing a. The /agn dependency on the host 
galaxy morphological type is clearly seen when early-type 
and late-type galaxies are divided. In general, /agn of 
early-type galaxies is much lower than that of late-type 
galaxies, especially at high luminosity and high a. 

At fixed luminosity, /agn of early- type galaxies is a de- 
creasing function of a, indicating that galaxies with lower 
velocity dispersion {<t< 200 km s~^) are dominant AGN 
hosts among early-types. The decreasing /agn with in- 
creasing a suggests that early-type galaxies with more 
massive black holes (assuming the Mbh — cr relation) are 
less likely to host AGN. Presumably, these most massive 
galaxies and their black holes are already well evolved as 
dormant black holes in the present-day universe. In con- 
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Fig. 7. — Distributions of AGN host {blue dots) and all galaxies in the sample {gray dots) in the luminosity versus velocity dispersion 
plane. The contours denote constant AGN fractions (colored solid lines) and galaxy number densities {dotted lines) where the bin sizes are 
A{loga) = 0.02 and A(Mr)=0.14. The galaxy density contour levels are on logarithmic scales. 




u — r 



Fig. 8. — The fraction of AGNs in early-type (solid lines) and late-type galaxies (dotted lines) as a function of luminosity {upper) and 
velocity dispersion (lower). The sizes are A(Mr)=0.14, A{logo-) = 0.02, and A(u — r) = 0.08. 
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trast, late-type galaxies with intermediate a (~130 km 
s~^) show highest /agn- This is because the fraction of 
red galaxies increases as a increases, and that /agn is 
actually a monotonically increasing function of a in the 
case of late- type galaxies (see §4.2). 

At fixed cr, /agn increases with luminosity for early- 
type galaxies, while this dependency becomes weaker for 
late-type galaxies. The decrease of /agn at lower lumi- 
nosity and lower a for late-type galaxies is not straight- 
forward to interpret, however, we speculate that it could 
be due to the lack of black holes in some fraction of 
low-mass disk-dominated galaxies or the relatively weak 
AGN activity (e.g., low [OIII] luminosity for given Ed- 
dington ratios) in these galaxies. 

To investigate whether this trend is caused by some 
selection effect due to the [OIII] flux limit in the survey 
as discussed in §3.3, we used a subsample of high-power 
AGNs with L[oiii]/Mbh> 1 (Lq/Mq) in the calculation 
of the /agn- The fraction of the high-power AGNs is 
much lower than that of all AGNs. However, the overall 
trend is qualitatively the same and is not due to the [OIII] 
flux limit of the sample. 

4.2. Color 

In this section, we investigate how /agn is related with 
the host galaxy u — r color, which is an indicator of star 
formation history or mean stellar age. Again, we ex- 
clude galaxies with cr < 90 km s~^ to reduce systematic 
errors induced by the instrumental resolution of SDSS 
spectroscopy. However, the results remain qualitatively 
the same even when the whole sample is used. 

The /agn dependency on host galaxy color is presented 
in Mr or a versus the u — r color plane (Figure [8]) . In 
general AGN host galaxies show a narrower distribu- 
tion in M — r color than non-AGN galaxies, indicating 
that galaxies with intermediate color are dominant AGN 
hosts. The number density of AGN host galaxies peaks 
at u — r ~ 2.4. The high fraction of AGN host galax- 
ies in the intermediate-color region between the red se- 
quence and the blue cloud is consistent with other studies 
(e.g. Nandra et al. 2007; Schawinski et al. 2007; Mar- 
tin et al. 2007). For early- type galaxies, the peak of 
/agn is clearly offset from the peak of number density 
of non-AGN galaxies. The /agn of early- type galaxies 
monotonically increases as the galaxies become bluer and 
reaches the maximum at w — r ~ 2. Similarly, for late- 
type galaxies /agn is maximum at intermediate color 
u - r = 2 - 2.4. 

In contrast to the dependency on color, /agn of early- 
type galaxies does not strongly depend on luminosity 
or velocity dispersion. The middle column of Figure [8] 
shows that /agn of early-type galaxies depends almost 
entirely on color, implying that triggering of AGN activ- 
ity is independent of the luminosity or velocity dispersion 
at fixed galaxy color. The dependence of /agn on Mr 
or (7 found in the previous section is merely due to the 
fact that Mr or a is correlated with m — r in the sense 
that early-type galaxies with higher velocity dispersion 
are redder and that for given tr, luminosity is higher for 
AGN host galaxies than non-AGN galaxies (see Figure 4 
upper panel). 

Combined with Figure [7l these results indicate that 
/agn strongly depends on host morphology and u — r 
color. For early-type galaxies, it does not much depend 



on luminosity or velocity dispersion. Thus, one of the 
necessary conditions for early-type galaxies to trigger 
AGN activity is the bluer color, which is an indication of 
the presence of young stellar population and/or cold gas 
that can be accreted to the black hole. In contrast, late- 
type galaxies show weak net dependence of /agn on Mr 
and a in addition to u — r (right panels in Fig. [5]) . While 
late-type galaxies with intermediate color show the high- 
est /agn, /agn at fixed color increases with luminosity 
or velocity dispersion. The dependence on luminosity 
or velocity dispersion becomes much weaker for the red- 
der late-type galaxies (it — r > 2.6), as similarly shown 
among early-type galaxies. Note that for late-type AGN 
host galaxies, the range of color distribution becomes 
larger with increasing velocity dispersion, implying that 
at lower mass scales AGN activity is more difficult to 
be triggered and only occurs in narrower galaxy color 
ranges. 

4.3. Color gradient 

We investigate how /agn depends on the color gradient 
in FigurelH where the distribution of /agn is presented in 
the u—r color versus /S.{g — i) plane. Note that early-type 
and late-type galaxies occupy distinct regions because 
galaxy morphology is mainly determined in this plane 
(see Figure [3]). 

For early-type galaxies, the dependence of /agn on 
A((7 — i) is not strong. At u — r > 2.4, where most of 
early-type galaxies are located, the /agn suddenly drops, 
indicating that typical early-type galaxies are not likely 
to host AGN. At bluer colors {u — r < 2.4) where the 
/agn is highest, the /agn is independent oil^{g—i). The 
galaxies in this color range have on average more positive 
A(5 — i) (i.e. bluer center). Perhaps, this may indicate 
the association between AGN and blue-core spheroids 
(e.g. Menanteau et al. 2005) and imply that AGN ac- 
tivity is strongly correlated with the amount of cold gas 
or younger stellar population in the bulge. This trend 
is also found when the subsample of high-power AGNs ( 
-^[oiii]A^^BH> 1 Lq/Mq) is used. 

In the case of late-type galaxies, the AGN fraction has 
negligible dependence on A (17 — i) at fixed color, probably 
due to the overall presence of cold gas and star-formation 
among late-type galaxies. 

4.4. Light concentration 

We investigate how /agn varies with light concentra- 
tion (cin) in Figure[T0l The dependence of /agn on Cin is 
relatively weak compared to that of u—r color. For early- 
type galaxies, /agn shows only negligible dependence on 

Cin- 

At u — r ~ 2, the /agn becomes dependent on Qn, es- 
pecially for late-type galaxies. Among late-type galaxies 
with u — r < 2.6, more concentrated ones have higher 
/agn at a given color. On the other hand, for redder 
late type galaxies ( u — r > 2.6), /agn is independent 
of the light concentrations. The dependence of /agn 
on Cin is very similar to that on a as shown in Fig- 
ure |8] due to the correlation between Cin and a. Given 
the tight correlations between Cin and cr, and between 
Mbh and host galaxy bulge mass, these results may sug- 
gest that among late-type galaxies with enough cold gas 
(u — r < 2.6), more bulge-dominated galaxies hosting 
more massive black hole are more likely to trigger AGN 




u - r 



Fig. 9. — Distributions of AGN host (blue dots) and all galaxies {gray dots) in the u — r color versus g — i color gradient plane. The bin 
sizes are A(A(3 - i)) = 0.026 and A(n - r) = 0.08. 




u - r 

Fig. 10. — Distributions of AGN host [blue dots) and all galaxies {gray dots) in the concentration index versus color plane. The bin sizes 
are A(A(g - i)) = 0.01 and A{u - r) = 0.08. 



activity. For the redder late-types, which are on aver- 
age more bulge-dominated, the amount of cold gas seems 
more requisite ingredient for AGN activity. The trend re- 
mains the same when we use a subsample with a narrow 
(T range (100 to 170 km s^^) in order to separate the 
effect of Cin from that of a. 

5. EFFECTS ON AGN POWER 

In this section, we examine how the host galaxy prop- 
erties are related with the power of AGN activity, by 
adopting L[oiii] normalized by Mbh as an Eddington ra- 
tio indicator (see Fig. [5] and § I3.3p . First, we compare 
host morphology with AGN power (§ IS.ip . Then, we 
investigate how color, color gradient, and concentration 
index are related with AGN power (ij 15. 2p . 

5.1. Host morphology 

Figure [TT] presents AGN number density as a func- 
tion of AGN power (Ljqjjjj/Mbh)- Here, we consider 



AGNs with log Mbh/Mq > 6.7 (left panels) since the 
systematic error due to the instrumental resolution of 
SDSS spectroscopy becomes severe at (T<90km s'-^. Re- 
gardless of host morphology, the Eddington ratio ranges 
over 3 orders of magnitude. The distribution of AGNs 
hosted by late-type galaxies peaks at ^[oiiij/Mbh ^ 0.3 
(Lq/Mq). For AGNs hosted by early- type galaxies, the 
peak of the distribution is slightly lower by ~0.2 dex. 
The large range of the Eddington ratio distribution of 
Type II AGNs is consistent with that of broad-line AGNs 
(Woo & Urry 2002; cf. KoUmeier et al. 2006), indicating 
that AGN power ranges from the Eddington luminosity 
to 10~^ of the Eddington limit. 

The number density decreases at the high end of the 
distribution, reflecting an intrinsic luminosity function at 
fixed Mbh- In contrast, the decrease of the number den- 
sity at the low end of the distribution is probably due to 
the selection effects since weak [OIII] lines are hard to de- 
tect for a given measurement limit of the [OIII] line flux. 
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Fig. 11. — The Eddington ratio distributions (upper) and rel- 
ative fraction of host galaxy morphology (lower) for AGNs with 
log Mbh/Mq > 6.7(left) and AGNs with 7.5 < log A/bh A-f < 8 
(right). 



The number density may keep increasing down to non- 
AGN level if we include objects with very low level AGN 
activity, for which emissions from HII regions in the host 
galaxy dominate over AGN emission lines and/or [QUI] 
line is intrinsically too weak to detect with a few A EW 
detection limit. This line of conjecture can be supported 
by the fact that the overall AGN fraction in our sample is 
less than 10%. This is much lower than the result of the 
nearby galaxy survey by Ho et al. (1997), where many 
dwarf AGNs with weak emission lines (EW-0.25A) were 
detected owing to the high spatial resolution ^200 pc in 
their survey. Thus, the AGN number density in our study 
should be taken as a lower limit since our AGN sample do 
not include very weak emission line AGNs. The intrinsic 
distribution of the Eddington ratio is closely related with 
the duty-cycle of AGN activity and beyond the scope of 
this paper. 

The striking difference between AGNs hosted by galax- 
ies of different morphology comes from the overall shape 
of the Eddington ratio distributions. Late-type galaxies 
are dominant hosts over all AGN power, while the frac- 
tion of early-type host galaxies increases for low-power 
AGNs (L[o„i]/AfBH < 0.1, or < 1% of the Eddington 
limit). 

As shown in Figure [5l the Eddington ratio range of 
lower Mbh objects is biased toward high-power due to 
the detection limit of [QUI] line (i.e. log Ljomj/i© 
^ 5 — 6). Hence, the peak of the distribution observed 
in the total sample may have shifted toward higher Ed- 
dington ratios. To avoid this selection effect and secure 
a consistent range of the Eddington ratio, we consider 
AGNs in a narrow Mbh range, 7.5 < log Mbk/Mq < 
8 (right panels in Fig. [TT|l . This narrow range of Mbh 
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Fig. 12. — Galaxy properties: u — r color (upper), color gradi- 
ent (middle), and concentration index (lower) as a function of the 
Eddington ratio. AGNs with 7 < logMBn/M© < 8 are plotted to 
avoid selection effects. The error bars include measurement errors 
of each parameter as well as the scatter in each bin. 

is also useful to avoid any systematic effect caused by 
slightly different Mbh distributions of two morphological 
subsamples. Figure [11] right panels show that the Ed- 
dington ratio distributions of AGNs in the narrow Mbh 
ranges show the same trend, confirming that early-type 
galaxies host on average less powerful AGNs and that the 
relative fraction between early-type and late-type hosts 
is a strong function of AGN power. 

5.2. Color, color gradient, and concentration 

Figure [12] presents host galaxy color, color gradient, 
and concentration as a function of AGN power. AGNs 
with 7 < logMBn/Af© < 8 are presented here to mini- 
mize selection effects since for lower mass AGNs the Ed- 
dington ratio distribution is biased against low-power ob- 
jects (see Fig.[S]). In this Mbh range, the Eddington ratio 
distribution is not heavily biased down to i[oiii]/-^-^BH 
- 0.1 (~1% of the Eddington hmit), and the [OIII] lu- 
minosity is typically larger than 10^ Lq. 

Top panels in Figure [T^] present the color distribution 
of host galaxies as a function of AGN power, showing that 
low-power AGNs are hosted by slightly redder galaxies 
especially in late-type galaxies, although at given AGN 
power there is a spread color distribution. As AGN power 
increases from 1% to 100% of the Eddington limit, the 
average u — r color of late-type galaxies becomes bluer 
by ~ 0.2. 

In the case of color gradient (middle panels), AGN 
power do not show strong dependence since galaxy colors 
are averaged in the color gradient space. A large frac- 
tion of blue AGN hosts tend to have bluer inner part (as 
shown in Fig. [9]), and have higher Eddington ratios, how- 
ever, this trend seems insignificant in the color gradient 
space. 

The concentration index also shows no strong trend. 
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Among early-type hosts, slightly less concentrated galax- 
ies host more powerful AGN. However, as in the case of 
color gradient, the range in the concentration index is 
large at any given AGN power, indicating that this de- 
pendency is weak. As discussed in Section 4, when the 
colors of galaxies are limited to a small range so that c-m 
or A(g — i) becomes a determining factor of /agn, its de- 
pendence on AGN power is slightly improved, however, 
it is still not significant. 

6. SUMMARY AND DISCUSSION 

We conservatively identified 11,521 Type II AGNs in a 
large volume- and luminosity-limited sample of 144,940 
galaxies with 0.025 < z < 0.107, selected from SDSS 
DR5, to investigate how host galaxy properties arc re- 
lated with AGN fraction and power. The main results 
are summarized as follows: 

1. Among other galaxy properties, u — r color shows 
the dominant dependence of the AGN fraction. In 
particular for early-type galaxies, the /agn does 
not depend on luminosity or velocity dispersion at 
fixed u — r color. 

2. AGNs arc typically hosted by intermediate-color 
late-type (m — r = 2 ~ 2.4) and bluish early-type 
galaxies (peak at u — r ~2.0), indicating that AGN 
host galaxies have younger stellar population than 
non-AGN galaxies at given luminosity or velocity 
dispersion. 

3. AGNs arc dominantly hosted by latc-type galax- 
ies with intermediate velocity dispersion (~130 km 
s~^), while among galaxies with highest velocity 
dispersion {a > 200 km s~^), the number of early- 
type host galaxies is comparable to that of late- type 
galaxies. 

4. Regardless of morphological types, the fraction of 
AGN increases with galaxy luminosity when their 
velocity dispersions are fixed. In contrast, at fixed 
luminosity the /agn of early-type galaxies mono- 
tonically decreases as velocity dispersion increases, 
while the /agn of late-type galaxies peaks at inter- 
mediate velocity dispersions {a 130 km s~^). 

5. The AGN fraction is independent of the color gra- 
dient in the case of late- type and bluer (u — r < 2.4) 
early-type galaxies. 

6. Among late-type galaxies, more bulge-dominated 
(higher velocity dispersion and more concentrated) 
galaxies show a higher AGN fraction. 

7. Late-type galaxies are dominant host of AGN of all 
power (10~^ < Eddington ratio) while the fraction 

of early-type host galaxies increases at low power 
(< 1% of the Eddington limit). For both mor- 
phological types, the Eddington ratio derived from 
-^[Oiii]/-^BH ranges over 3 orders of magnitude, in- 
dicating various levels of accretion for given Mbh- 

8. On average bluer late- type galaxies host more pow- 
erful AGNs while color gradient and concentra- 
tion index do not show strong dependence on AGN 
power. 



Our results indicate that the requisite ingredient for 
triggering AGN activity is massive (high velocity dis- 
persion) bulges and intermediate colors. The necessary 
condition for AGN activity seems mostly different with 
galaxy morphology. These results are consistent with 
a scenario that in the present-day universe, early-type 
galaxies with higher velocity dispersion and redder color 
are harder to host AGNs since these galaxies already con- 
sumed gas at the center or do not have sufficient gas sup- 
ply to the central black hole. In contrast, intermediate- 
velocity dispersion, intermediate-color, and more concen- 
trated late- type galaxies are more likely to host AGNs. 
Perhaps, this trend results from that some fraction of 
low-mass, blue, and less concentrated late-type galaxies 
may not host massive black holes, as predicted by black 
hole seed models (e.g. Volontary et al. 2008) or may 
host very low-power AGNs. The presence of black holes, 
either dormant or active, in low-mass disk-dominated 
galaxies is currently a topic of active research (see re- 
lated discussions by Gallo et al. 2008; Ferrarese et al. 
2006) 

It is not straightforward to compare our results with 
previous studies due to several differences in the sample 
selection and analysis. The different schemes we adopted 
were to match galaxy properties when we compared AGN 
host and non-AGN galaxies, and to normalize the AGN 
luminosity by black hole mass. While Kauffmann et al. 
(2003) found that AGNs reside almost exclusively in mas- 
sive galaxies, our results show that late-type galaxies 
with intermediate-luminosity and intermediate-velocity 
dispersion are dominant AGN host in the present-day 
universe. This different result is probably due to the fact 
that our AGN sample is more conservatively defined with 
a higher S/N cut {S/N > 6) on each emission lines used 
in flux ratio diagrams, and do not include very low-power 
AGNs, which are perhaps more prevalent among massive 
galaxies. 

Although we find strong dependence of AGN activity 
on host galaxy properties, we do not find any sufficient 
host galaxy property to trigger AGN activity since there 
are more non-AGN galaxies than AGN host galaxies for 
any given combination of host galaxy properties. These 
findings may imply that large scale environments are cru- 
cial for triggering AGN activity or that local processes at 
the vicinity of the central black holes play an important 
role, especially for the dominant low-luminosity AGNs in 
the present-day universe. It is also possible that the life 
time of AGN is shorter than that of imprint of triggering 
mechanisms so that distinct galaxy properties, i.e.. inter- 
mediate color, are still observable although AGN activity 
finished. We note that these findings are relevant only 
for local galaxies and that host galaxy properties may 
show very different relations with AGN activity at high 
redshifts. 
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